REPORT  DOCUMENTATION  PAGE 

_ .  .  afrl-sr-bl-tr-01 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services 

Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperv\ 

IS, 

iis 

ozh  2,  __I 

1.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPOR.  . . 

Sep  4,  2000  Final  Report  1 

„  vnlbU ouvcncu 

12  April  2001-30  November  2000 

4.  TITLE  AND  SUBTITLE 

Global/Local  Design  Optimization  of  A  Power  Distribution  System 

5.  FUNDING  NUMBERS 

F49620-00- 1-0056 

6.  AUTHOR(S) 

Zafer  Gurdal,  Scott  Ragon 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Virginia  Polytechnic  Institute  &  State  University 

301  Burruss  Hall 

Blacksburg,  VA  24061-0249 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

AFOSR/NM 

801  N.  Randolph  St,  Rm  732 

Arlington,  VA  22203-1977 

10.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER 

F49620-00- 1-0056 

11.  SUPPLEMENTARY  NOTES 


12a.  DISTRIBUTION  AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 


AiR  K))!CE  OFFICE  OF SeMfr  8B8BWIGtttWi9BP>;'opE~ 
NOTICE  OF  TRANSWITTAL  Q 
HAS  BEEN  REVIEWED  AND 
LAWAFR 180-12.  DISTRIBliili 


TiC.  THIS  TECHNICAL  REPORT 
S  APPROVED  FOR  PUBUC  RELEASE 
ION  IS  UNLIMITED. 


13.  ABSTRACT  (Maximum  200  words) 

The  MATLAB/Simulink  analysis  tools  developed  and  upgraded  by  the  electrical  engineering  team  were  successfully 
incorporated  into  the  optimization  environment.  A  well  posed  optimization  problem  was  formulated  for  a  sample  system 
consisting  of  a  filter  and  converter.  It  was  shown  that  the  solution  of  the  optimzation  problem  yielded  realistic  parameter 
values.  Also  the  existence  of  suboptimal  solutions  was  investigated.  These  results  are  the  core  results  for  developing  an 
optimization  methodology  for  the  full  power  distribution  system.  The  optimization  methodology  is  being  reworked  into  a 
global/local  formulation.  This  formulation  allows  us  to  define  an  optimization  problem  for  the  entire  power  distribution 
system.  A  preliminary  framework  for  the  sample  system  has  been  developed,  and  the  appropriate  software  is  being 
generated.  The  global/local  formulation  has  also  been  applied  to  the  design  filters  components  with  the  inductors  as  their 
local  level  subcomponents.  Preliminary  results  indicate  a  significant  decrease  in  computational  time  while  retaining  the 


accuracy  of  the  full  problem. 


20010501  0/5 


14.  SUBJECT  TERMS 


15.  NUMBER  OF  PAGES 

13 


16.  PRICE  CODE 


17.  SECURITY  CLASSIFICATION 
OF  REPORT 


18.  SECURITY  CLASSIFICATION 
OF  THIS  PAGE 


19.  SECURITY  CLASSIFICATION 
OF  ABSTRACT 


Standard  Form  298  (Rev.  2-89}  (EG) 

Prescribed  by  ANSI  Std.  239.18 

Designed  using  Perform  Pro,  WHS/DIOR,  Oct  94 


UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


nM  APR  1  2  2001 


Global/Local  Design  Optimization  of  A  Power  Distribution  System 


Final  Report 


Sponsored  By 

Air  Force  Office  of  Scientific  Research 
AFOSR  Contract  Number:  F49620-00-1-0056 

Project  ending  date:  1 1/30/00 


Dr.  Neal  D.  Glassman,  Program  Manager 
AFOSR/NM 

801  North  Randolph  Street,  Room  732 
Arlington,  VA  22203-1977 


September  4,  2000 

Prepared  by 

Zafer  Gurdal,  Scott  Ragon 

Engineering  Science  and  Mechanics  Department  -  0219 
in  cooperation  with 


Douglas  Lindner 

Bradley  Department  of  Electrical  and  Computer  Engineering  -0111 

Virginia  Tech 
Blacksburg,  VA  24061 


Objectives 


Assist  in  the  development  of  optimization  methodologies  for  the  design  of  next  generation  power  systems. 
Implement  a  global/local  optimization  technique  to  achieve  a  cost  effective  design  optimization  strategy. 

Status  of  the  Effort 

During  the  past  year  significant  progress  was  made  on  demonstration  of  the  global/local  methodology  for 
minimum  weight  design  of  the  power  electronic  circuits. 

Accomplishments/New  Findings 

1.  The  MATLAB/Simulink  analysis  tools  developed  and  upgraded  by  the  electrical  engineering  team  were 
successfully  incorporated  into  the  optimization  environment. 

2.  A  well-posed  optimization  problem  was  formulated  for  a  sample  system  consisting  of  a  filter  and  a  converter.  It 
was  shown  that  the  solution  of  the  optimization  problem  yielded  realistic  parameter  values.  Also  the  existence  of 
suboptimal  solutions  was  investigated.  These  results  are  the  core  results  for  developing  an  optimization 
methodology  for  the  full  power  distribution  system. 

3.  The  optimization  methodology  is  being  reworked  into  a  global/local  formulation.  This  formulation  allows  us  to 
define  an  optimization  problem  for  the  entire  power  distribution  system.  A  preliminary  framework  for  the  sample 
system  has  been  developed,  and  the  appropriate  software  is  being  generated.  The  global/local  formulation  has  also 
been  applied  to  the  design  filters  components  with  the  inductors  as  their  local  level  subcomponents.  Preliminary 
results  indicate  a  significant  decrease  in  computational  time  while  retaining  the  accuracy  of  the  full  problem. 
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Publications 

The  following  journal  article  was  submitted  for  publication. 

Chandrasekaran,  S.,  S.  A.  Ragon,  D.  K.  Lindner,  Z.  Gurdal,  and  D.  Boroyevich,  "Optimization  of  a  Power 
Distribution  Subsystem,"  submitted  to  AIAA  Journal  of  Aircraft. 

The  following  conference  paper  has  been  presented  (name  of  the  presenter  is  highlighted). 

Ragon,  S.  A.,  S.  Chandrasekaran,  Z.  Gurdal,  D.  K.  Lindner,  and  D.  Boroyevich,  "Optimal  Design  of  a  Power 
Distribution  System,"  accepted  for  the  38th  AIAA  Aerospace  Sciences  Meeting  and  Exhibit:  Multidisciplinary 
Design  Optimization,  Reno,  NV,  Jan  10  -  13,  2000. 


Interactions/Transitions 


7a.  Participation  at  meetings,  etc 

38th  AIAA  Aerospace  Sciences  Meeting  and  Exhibit:  Multidisciplinary  Design  Optimization,  Reno,  NV,  Jan  10  - 
13, 2000:  Prof  Giirdal  presented  a  paper  on  optimization  of  the  sample  system. 

Exergy  Workshop,  AFRL,  Dayton,  OH,  Dec  16,  17,  1999:  This  workshop  was  organized  by  AFRL  to  explore  the 
possibility  of  using  exergy  as  a  basis  for  the  optimization  of  various  aircraft  subsystems.  Profs  Giirdal  and  Lindner 
presented  the  results  on  optimization. 

Exergy  Workshop,  Virginia  Tech,  Blacksburg,  VA,  March  23-24,  2000:  This  workshop  was  a  follow-up  to  the 
previous  Exergy  Workshop.  Prof.  Lindner  presented  results  on  the  global/local  formulation  for  a  power 
distribution  system. 

7c.  Transitions 

Based  on  the  preliminary  results  for  the  optimization  of  the  sample  system,  the  Center  for  Power  Electronic 
Systems  (CPES)  at  Virginia  Tech  was  awarded  a  contract  by  Schnieder  Electric  Company  to  (in  part)  develop  an 
optimization  tool  for  the  design  of  a  power  factor  correction  circuit.  (This  circuit  is  basically  a  filter  followed  by  a 
boost  converter  as  in  the  sample  system  studied  under  this  project.)  Profs.  Boroyevich,  Giirdal,  and  Lindner  are 
among  the  Pi’s  for  this  project.  Mr.  M.  Arpilliere  of  Schnieder  Electric,  S.A.  is  the  Program  Manager.  The 
optimization  methodology  developed  under  the  AFOSR  project  was  extended  to  the  power  correction  circuit. 

Additionally,  under  the  same  contract,  ADOPTECH  Inc,  of  Blacksburg  Virginia  initiated  a  preliminary  study  for 
implementation  of  the  Genetic  Algorithms  for  power  electronic  circuit  design.  A  commercial  code  capable  of 
handling  discrete  power  electronics  problems  with  specified  topology  is  being  developed.  An  expansion  of  the 
capability  to  design  power  electronics  circuit  topologies  is  envisioned  if  the  theoretical  developments  from  the 
Virginia  Tech  team  can  be  realized. 


New  Discoveries,  Inventions,  or  Patent  Disclosures 

There  were  no  new  disclosures  during  this  period. 

Honors/Awards 

No  honors  or  awards  were  accepted  during  this  period. 
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Introduction 


Global/Local  Design  Optimization 


A  major  goal  of  the  present  research  project  is  to  develop  mathematical  optimization  methodologies  and  tools  for 
the  design  of  next  generation  power  systems  for  aircraft/spacecraft.  The  design  task  is  quite  complex,  as  power 
distribution  systems  consist  of  many  diverse  subsystems,  each  of  which  are  typically  governed  by  complex  sets  of 
non-linear  equations.  The  simulation  of  the  global  power  distribution  system  involves  complicated  interactions 
between  the  various  subsystems  that  need  to  be  accounted  for  in  the  design  process.  Our  earlier  work  has 
demonstrated  that,  even  for  a  very  simple  two-component  circuit  design,  it  is  advantageous  to  optimize  the  system 
as  a  whole  rather  then  optimizing  each  of  the  individual  components  independently. 

As  more  and  more  components  are  added  to  the  optimization  problem,  certain  difficulties  arise.  The  optimization 
algorithm  becomes  less  efficient  as  the  number  of  design  variables  and  constraints  increase,  and  the  computational 
expense  of  simulating  the  response  of  the  overall  power  system  rapidly  grows  very  large.  In  order  to  surmount 
these  difficulties,  we  have  formulated  a  global/local  design  methodology  for  the  power  system  design. 

In  the  global/local  methodology,  the  overall  design  problem  is  decomposed  into  two  design  levels:  a  global  design 
level  and  a  local  design  level.  At  the  global  design  level,  the  overall  system  is  designed  using  a  set  of  global  design 
variables.  These  global  variables  do  not  specify  the  details  of  each  component,  but  instead  describe  the  effect  of 
each  component  on  the  other  components  in  the  system.  The  local  design  level  consists  of  a  number  of  separate 
detailed  design  models  of  each  individual  component  type.  For  example,  there  may  be  a  detailed  design  model  of 
an  input  filter,  a  detailed  design  model  of  a  DC-DC  buck  converter,  etc.  The  local  design  variables  specify  the 
details  that  are  neglected  at  the  global  level. 

The  global  and  local  design  levels  interact  with  one  another  through  a  local  design  database.  Before  the  system 
level  design  process  begins,  a  large  number  of  optimized  local  designs  of  each  power  system  component  are 
generated  and  stored  in  a  design  database.  Each  optimized  local  design  is  generated  as  a  function  of  a  set  of  global 
parameters  such  as  terminal  voltages/currents  and  input/output  impedances.  Once  the  design  database  has  been 
constructed,  it  is  used  by  the  global  optimizer  to  design  the  overall  power  system.  At  each  global  design  iteration, 
the  global  optimizer  is  able  to  consult  the  design  database  to  determine  the  weight  of  an  optimized  component  as  a 
function  of  the  globally  determined  parameters.  In  this  way,  the  global  optimizer  can  incorporate  information  about 
optimized  local  designs  without  actually  performing  a  local  optimization  at  each  iteration. 

A  more  detailed  description  of  the  design  formulation  of  a  global/local  power  distribution  system  will  be  presented 
below.  First,  however,  the  basic  global/local  design  methodology  will  be  demonstrated  using  a  simpler  power 
electronics  design  problem  that  we  refer  to  as  the  “global/local  inductor  design  problem”. 


Global/Local  Inductor  Design  Problem 

In  this  section,  the  global/local  design  problem  will  be  demonstrated  not  for  the  system  level  design  problem,  but 
for  a  simpler  power  electronics  design  application.  This  simpler  application  will  enable  us  to  evaluate  the  process 
without  loosing  the  generality  of  it  when  applied  to  the  overall  system  design.  This  application  is  concerned  with 
the  design  of  the  input  filter  illustrated  in  Figure  1:  the  three  inductors  Lj,  L2 ,  and  Z^are  designed  in  detail  at  the 
local  level,  and  the  filter  as  a  whole  is  designed  at  the  global  level. 

In  our  previous  work,  all  of  the  design  variables  describing  the  filter  (see  Table  1)  were  considered  simultaneously. 
Note  that  a  large  number  of  these  variables  (15  variables  out  of  the  18  in  the  table)  describe  physical  dimensions 
and  characteristics  of  the  three  inductors.  For  example,  the  variables  ndi  Acpd,  C^,  Wwd,  and  lgd  all  specify  the 
detailed  design  of  inductor  Ld. 


Number  of  turns  for  L2 


Cross  sectional  area  of  winding  for  L2 


Center  leg  width  for  L2 


Window  width  for  L2 


Airgap  length  for  L2 _ 


If  the  filter  design  problem  is  examined  in  more  detail,  it  can  be  seen  that  the  problem  naturally  lends  itself  to 
decomposition  into  global  and  local  design  levels.  If  we  wish  to  determine  the  behavior  of  the  filter  as  a  whole,  it  is 
not  necessary  to  know  the  details  of  each  inductor  design:  we  only  need  to  know  the  inductance,  L ,  of  each 
inductor.  This  indicates  that  we  may  be  able  to  remove  the  design  variables  describing  the  inductor  details  from  the 
overall  (global)  design  problem.  If  we  wish  to  compute  the  total  weight  of  the  filter,  however,  we  need  to  know  the 
weight  of  each  inductor  as  a  function  of  L.  As  it  turns  out,  it  is  possible  to  formulate  a  separate  local  design 
problem  that  will  allow  us  to  find  the  minimum  inductor  weight  as  a  function  of  the  inductance,  L,  and  the  peak 
inductor  current,  IL(pk).  This  second  optimal  design  problem  will  constitute  the  local  design  level  for  the  filter 
design  problem,  and  can  be  stated  as  follows: 


Local  Inductor  Design  Problem 


For  specified  parameters:  L,  IL(Pk) 
Minimize:  inductor  weight 
Using  variables:  «,  Acp,  Cw,  Ww,  lg 


Such  that: 

•  The  computed  inductance: 

r 

i  +Z  /fi’ 

X  Pi 

is  equal  to  the  given  inductance,  L.  Here,  fX0  is  the  permeability  of  free  space,  pr  is  the  permeability  of 
ferrite,  K  <  is  the  aspect  ratio  of  the  core  leg,  and  Zp  is  the  mean  magnetic  path  length. 

•  The  widths  of  the  center  leg,  Cw,  and  of  the  window,  Ww,  are  greater  than  1  mm. 

•  The  copper  wire  cannot  have  a  cross-sectional  area  less  than  7.29  x  10"8  m2. 

•  The  current  density  in  the  windings  of  the  inductor  cannot  be  greater  than  maximum 
allowable  current  density  for  copper  which  was  considered  to  be  5.5  x  106  /1/m2. 

•  The  window  of  the  EE-core  houses  the  windings  and  the  bobbin.  The  area  occupied  by 
the  windings  is  given  by: 


nA 

W  = — — 
a  F  ’ 

W 

where  Fw=0A  is  the  window  fill  factor.  The  window  fill  factor  is  included  to  account  for  imperfections 
in  the  windings  on  the  bobbin.  The  area  occupied  by  the  bobbin  in  the  window,  using  simple 
geometry,  can  be  determined  as: 

Wb=WbobK2Ww, 

where  Wbob  is  the  thickness  of  the  bobbin  wall.  Hence,  the  window  is  required  to  be  large  enough  to 
accommodate  the  windings  and  the  bobbin.  This  requirement  is  formulated  as  a  constraint  given  by: 

K2Wl>Wa+Wb. 


•  In  order  to  prevent  the  inductor  core  from  running  into  saturation,  the  dimensions  of 
the  inductor  should  be  such  that  the  maximum  allowable  saturation  flux  density  for  a 
ferrite  core  material,  Bsp  =  0.3  T,  is  not  exceeded.  The  maximum  flux  density  is 
determined  as  the  ratio  of  the  maximum  flux  to  the  area  of  cross  section  of  the  center 
leg.  Hence,  this  constraint  is  given  by: 

d  .  *pt  _  (vW”)_  iLIL(pk)/n) 

'  Ac*  KXC2W  KXC2W 


By  solving  the  local  inductor  design  problem,  we  can  compute  the  minimum  inductor  weight  as  a  function  of  L  and 
IL(p k).  We  can  now  formulate  a  design  problem  for  the  filter  as  a  whole  in  which  we  replace  the  detailed  inductor 
design  variables  with  the  variables  Lu  L2 ,  and  Ld .  This  global  level  design  problem  can  be  stated  as  follows: 


Global  Filter  Design  Problem 
Minimize:  total  filter  weight 


Using  variables:  L,,  L2,  Ld,  C,,  C2,  Rd 

Such  that: 

The  magnitude  of  the  input-output  transfer  function  of  the  filter,  Av(jco)  remains 
between  an  upper  and  lower  bound  up  to  a  passband  edge  frequency,  6}  (see  Figure 

-1  dB<jAv(jO))f<6  dB  for  0  <co<cop  =2tt5x103 rad/sec. 

The  magnitude  of  the  input-output  transfer  function,  Av(jco),  remains  below  an  upper 
bound  above  a  certain  frequency,  ^  (see  Figure  2): 

KC/*>)|  <  ~60  dB  for  (o>cosX=2n- 50xl03 rad/sec 


Figure  2  -  Definition  of  frequency  domain  constraints 


For  a  given  load  disturbance,  the  maximum  transient  excursion  of  the  output  voltage 
of  the  filter,  AvB,  is  less  than  20  V. 

The  output  impedance  of  the  filter,  Zofi  which  depends  on  the  design  parameters,  must 
be  less  than  the  input  impedance  of  the  converter,  ZiL,  in  order  to  ensure  stability  of  the 
interconnected  system.  Since  ZiL  is  given,  sufficient  separation  between  these  two 
impedances  can  be  enforced  by  requiring  that  the  maximum  output  impedance  of  the 


=  max 


(M  b 


co 


<15  dBO  <Z(i. 


In  order  to  assure  that  the  interaction  between  the  filter  and  the  power  bus  is  stable 

'nput  Im^dance  of  the  filter  must  exceed  the  maximum  output 
mpedance  of  the  power  bus,  ZoC.  Since  ZoC  is  given,  sufficient  separation  between 

^f  the  mte^tisfy1'65  Can  bC  enf°rCed  ^  reqUiring  *at  the  minimum  mPut  impedance 
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Note  that  we  cannot  solve  the  global  design  problem  independently  from  the  local  design  problem  -  we  need  the 
local  design  problem  in  order  to  compute  the  weight  of  each  inductor  in  the  filter.  What,  then  has  been  gained  by 
decomposing  the  problem?  The  answer  lies  in  the  concept  of  the  local  design  database. 

Before  we  begin  solving  the  overall  filter  design  problem,  we  first  solve  the  local  optimization  problem  for  a  range 
of  different  L  and  IL(pk)  values,  and  store  the  results  in  a  local  design  database.  Given  the  data  stored  in  the 
database,  response  surface  methodologies  are  used  to  construct  a  simple  function  (usually  a  polynomial)  that 
approximates  the  inductor  weight  as  a  function  of  L  and  IL(pky 

Once  the  response  surface  function  is  obtained,  the  overall  filter  can  be  designed.  Whenever  the  weight  of  any  of 
the  three  inductors  is  required,  the  current  values  of  the  L  and  IL(pk)  are  inserted  into  the  response  surface  function 
and  the  required  weight  is  efficiently  computed  without  the  need  for  optimizing  the  local  component  for  that 
combination  of  the  L  and  IL(pk).  Given  the  existence  of  the  response  surface  function,  the  number  of  design  variables 
in  the  present  filter  design  problem  is  reduced  from  eighteen  to  six.  Likewise,  the  number  of  constraints  is  reduced 
from  twenty  (five  for  each  of  the  three  inductors  +  the  five  from  the  filter  as  a  whole)  to  five.  This  means  that 
individual  filter  optimization  problems  can  now  be  solved  more  efficiently  and  reliably.  Filter  design  studies 
(which  may  entail  solving  a  large  number  of  filter  optimization  problems,  each  corresponding  to  different 
operating  conditions,  etc.)  are  more  efficient  as  well,  as  the  same  local  database  can  be  reused  for  each 
optimization  problem.  In  fact,  the  local  design  database  can  potentially  be  reused  for  any  electrical  component  that 
contains  inductors. 

Results 

In  this  section,  preliminary  results  illustrating  the  global/local  inductor  design  methodology  are  presented.  First, 
the  filter  optimization  problem  was  solved  without  any  global/local  decomposition.  These  results  are  presented  in 
Table  2  (these  results  differ  from  those  generated  and  reported  earlier  in  this  project  because  i-  a  different  current 
density  for  copper  was  used,  ii-  a  slightly  more  accurate  expression  for  L  in  terms  of  the  inductor  design  variable 
was  used).  It  should  be  noted  that  it  took  a  number  of  different  optimization  “runs”  to  obtain  these  results  -  the 
optimization  problem  was  very  slow  to  converge  especially  with  respect  to  the  detailed  inductor  design  variables. 
The  detailed  inductor  design  variables  in  Table  2  correspond  to  the  final  values  of  Lj=25.6  pH ,  3.5  fiH,  and 

Ld=\1.9  pH. 


Next,  the  local  inductor  design  problem  was  solved  for  different  combinations  of  L  and  h(pk)  and  the  results  were 
stored  in  a  local  design  database.  A  third  order  polynomial  in  two  dimensions  was  fit  to  the  data;  this  polynomial 
approximates  the  optimal  inductor  weight  as  a  fimction  of  L  and  IL(pk).  This  polynomial  was  then  used  in 
conjunction  with  the  global  design  problem  to  generate  the  optimal  filter  design  shown  in  Table  3. 

These  results  agree  closely  with  those  presented  in  Table  2,  and  were  obtained  with  considerably  less 
computational  effort.  In  the  near  future  additional  global/local  inductor  designs  will  be  generated,  and  the 
differences  in  computational  effort  between  the  full  optimization  (no  decomposition)  and  the  global/local 
optimization  methodologies  will  be  documented. 
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Table  2  -  Optimal  filter  design  (no  decomposition) 


Variabl 

e 

Name 

Value 

Cj 

7.82  jiF 

c2 

31.96  pF 

Rd 

21.95  /xP 

rid 

2.47 

ACpd 

.940  x  10-6 
m2 

Cwd 

.287  x  10'2 
m2 

Wwd 

.521  x  10'2 
m2 

Igd 

.504  x  10'3 
m2 

Nj 

3.00 

Acpi 

.144  x  10'5 
m2 

cwl 

.388  x  10'2 
m2 

Wwl 

.679  x  10'2  m 

Igl 

.953  x  10'3 
m2 

n2 

8.24 

Acp2 

.108  x  10'5 
m2 

CW2 

.681  x  10‘2 
m2 

Ww2 

.941  x  10'2  m 

h2 

.200  x  10‘2 
m2 

Weight 

0.151  kg 

Table  3  -  Optimal  filter  design  obtained  using  global/local  methodology 


Variabl 

e 

Name 

Value 

Ci 

7.79  gF 

c2 

30.54  pF 

Rd 

22.61  jiF 

Lj 

26.15  gH 

l2 

277.44 

gH 

Ld 

18.36  fiH 

weight 

0.153  kg 
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Global/Local  Formulation  for  System  Level  Design 

The  formulation  for  the  system  level  global/local  design  problem,  which  is  one  of  the  primary  tasks  of  this  project, 
is  be  presented  in  this  section.  Here,  the  global  design  level  will  constitute  the  power  system  as  a  whole,  and  the 
local  design  level  will  include  detailed  local  models  of  the  individual  power  system  components,  such  as  the  filter 
component  described  above.  For  purposes  of  presenting  the  global/local  formulation,  we  will  concentrate  on  the 
simplified  global  system  illustrated  in  Figure  3. 


Figure  3  -  Block  diagram  of  sample  system 


Global  Level  Design 

The  global/local  formulation  presented  in  this  section  utilizes  the  model  illustrated  in  Figure  4  for  the  global  level 
design  problem.  In  this  model,  each  of  the  local  components  (in  this  case,  the  input  filter  and  the  buck  converter)  is 
represented  using  two-port  models.  In  this  representation,  each  component  in  the  power  system  is  characterized  by 
four  transfer  functions:  Au  Av,  Zit  and  Z0.  Depending  on  the  form  of  each  of  these  functions,  any  of  the  components 
in  the  power  system  may  be  modeled. 


Input  Filter  Buck  converter  — 


Figure  4  -  Global  representation  of  power  system 
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In  the  case  of  the  input  filter,  the  transfer  functions  are: 

A  =A  =1±V^ 
'F  if  A  (s) 

7  -r  lif/gk 
A  /  \ 


y  —  IT  _ V  j 

A(5)  =  l+— +  -4, 

1  j  a®,  < 


and  therefore  the  global  design  variables  consist  of  the  parameters  KoF,  KziF,  coaF ,  fii0F>  g^zf?  ^/f,  and  QF.  In  the 
case  of  the  buck  converter,  the  transfer  functions  are: 


ab=kb 


1 + s/(om 


A  _  p 

A(,) 
l  +  g/o>— 
“8  -  A(j) 

z  *7  A(5) 

"  F  1+j/®- 


Z.  =- 


A(s)  =  l  +  — - —  +  -^-. 

&«*. 

and  the  global  design  variables  are  the  parameters  AA>(.  ft},,/;,  ow,  ft}//;,  ft},/;.  and  The  two-port  global 

models  will  be  capable  of  accurately  predicting  the  low  frequency  behavior  of  the  power  system,  including  the 
interface  voltages  vB  and  v„.  Global  constraints  will  include  upper  bounds  on  these  voltages: 

Av  <  Av 

Avs  <  AvS  nux. 

In  addition,  the  stability  of  the  interconnected  system  will  be  assured  by  imposing  a  constraint  on  the  input 
impedance  of  the  converter  and  the  output  voltage  of  the  fdter: 


Local  Level  Design 

The  input  filter  local  model  (Figure  1)  and  the  buck  converter  local  model  (Figure  5)  are  the  same  as  those  used  in 
our  previous  optimization  work.  The  objective  of  the  local  level  design  problems  will  be  to  minimize  the 
component  weight.  The  local  design  variables  for  the  filter  were  listed  in  Table  1 ,  and  the  local  design  variables  for 
the  buck  converter  are  listed  in  Table  4. 
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Figure  5  -  Buck  converter  model 

Most  of  the  constraints  imposed  on  the  local  level  design  problems  are  identical  to  those  used  in  our  previous  work. 
Constraints  imposed  on  the  filter  at  the  local  level  include  frequency  domain  constraints  on  the  input-output 
transfer  function  (passband  and  stopband  constraints),  and  physical  constraints  governing  the  design  of  the 
inductors.  Local  buck  converter  constraints  include  frequency  domain  constraints  on  the  audiosusceptibility  of  the 
converter  and  the  crossover  frequency  of  the  voltage  loop  gain.  Constraints  will  be  imposed  to  guarantee  the 
internal  stability  of  the  converter,  and  the  physical  constraints  governing  the  inductor  design  will  be  enforced. 

In  addition  to  these  performance  and  stability  constraints,  we  will  require  that  the  local  level  filter  or  converter 
design  match  the  behavior  of  the  corresponding  global  level  design  up  to  a  specified  frequency.  This  will  be 
accomplished  by  matching  the  higher  order  local  level  transfer  function  with  the  lower  order  global  level  transfer 
function  as  closely  as  possible  up  to  this  frequency. 

Table  4  -  Design  variables  for  buck  converter 


Design  Variable 

Description 

N 

Number  of  turns  for  L 

Acp 

CSA  of  winding  for  L 

cw 

Center  leg  width  for  L 

Ww 

Window  width  for  L 

k 

Airgap  length  for  L 

c 

Capacitance 

fs 

Switching  frequency  of  buck  converter 

h0 

Voltage  controller  gain 

Voltage  controller  zero 

Voltage  controller  pole 

hip 

Current  controller  proportional  gain 

hu 

Current  controller  integral  gain 

For  the  system  level  design  problem,  a  different  local  design  database  will  be  constructed  for  each  different 
component  type  in  the  system.  In  the  present  case,  one  local  design  database  will  be  developed  for  the  filter,  and  a 
second  database  will  be  developed  for  the  converter.  Each  database  will  consist  of  number  of  optimized  local 
designs,  each  of  which  corresponds  to  different  combinations  of  the  global  design  variables.  The  local  filter 
database  will  be  constructed  as  a  function  of  the  global  filter  design  variables,  and  the  local  converter  database  will 
be  constructed  as  a  function  of  the  global  converter  design  variables. 
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Once  these  databases  are  constructed,  response  surface  methodologies  will  be  utilized  to  construct  functional 
approximations  to  the  minimum  weight  of  each  component  as  a  function  of  the  global  design  variables.  Once  these 
functions  are  developed,  it  will  be  possible  to  optimize  the  interconnected  system  in  an  efficient  manner.  For  any 
combination  of  the  global  design  variables,  these  functions  will  return  the  minimum  weight  of  the  local  component 
that  satisfies  all  local  constraints.  The  result  of  the  converged  optimization  process  will  be  a  minimum  weight 
system  that  satisfies  all  global  and  local  constraints.  We  will  continue  this  research  by  developing  local  design 
databases  for  both  the  filter  and  converter.  The  global/local  methodology  will  then  be  demonstrated  on  the  sample 
filter-buck  converter  system. 

In  addition,  the  local  level  design  optimization  problems  proved  to  be  mixed  discrete  continuous  optimization 
problems  with  predefined  set  of  gage  dimensions  (such  as  the  wire  dimensions),  off-the  shelve  component  selection 
(such  as  resistors  and  capacitors),  and  multiple  discrete  choices  (such  as  the  core  types  and  core  material  selection 
for  the  inductor  component,  and  different  component  types).  Such  problems  require  special  algorithms,  such  as  the 
genetic  algorithms,  to  be  adapted  to  the  power  circuit  design  optimization.  Finally,  some  of  the  local  component 
topologies,  such  as  the  filter  topology,  may  need  to  be  included  into  the  design  optimization  problem.  Such  an 
approach  will  again  require  a  powerful  search  scheme  capable  of  sifting  through  various  topologies  to  find  an 
optimal  solution.  An  extension  of  this  work  would  be  to  design  local  topologies  within  the  framework  of 
global/local  design  methodology  by  using  genetic  algorithms. 


